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ABSTRACT 



Context. The H ll region/PDR/molecular cloud complex S 106 is excited by a single 0-star The full extent of the warm and dense gas 
close to the star has not been mapped in spectrally resolved high-J CO or [C ll] lines, so the kinematics of the warm, partially ionized 
gas, are unknown. Whether the prominent dark lane bisecting the hourglass-shaped nebula is due solely to the shadow cast by a small 
disk around the exciting star or also to extinction in high column foreground gas was an open question until now. 
Aims. To disentangle the morphology and kinematics of warm neutral and ionized gas close to the star, study their relation to the bulk 
of the molecular gas, and to investigate the nature of the dark lane. 

Methods. We use the heterodyne receiver GREAT on board SOFIA to observe velocity resolved spectral lines of [C ll] and CO 1 1^10 
in comparison with so far unpublished submm continuum data at 350 yum (SHARC-II) and complementary molecular line data. 
Results. The high angular and spectral resolution observations show a very complex morphology and kinematics of the inner SI 06 
region, with many different components at different excitation conditions contributing to the observed emission. The [C ll] lines are 
found to be bright and very broad, tracing high velocity gas close to the interface of molecular cloud and Hll region. CO 11— >10 
emission is more confined, both spatially and in velocity, to the immediate surroundings of S106 IR showing the presence of warm, 
high density (clumpy) gas. Our high angular resolution submm continuum observations rule out the scenario where the dark lane 
separating the two lobes is due solely to the shadow cast by a small disk close to the star The lane is clearly seen also as warm, high 
column density gas at the boundary of the molecular cloud and H ll region. 
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1. Introduction 

The H II region S106 at a distance of 1.2 to 1.8 kpc 
dSchneider et al.L 12007) in Cygnus is a prominent bipolar emis- 
sion nebula associated with an extended molecular cloud. The 

■ nebula is excited by the single, late O-type star SI 06 IR, which 
creates a bright Photon Dominated Region (PDR) at the molecu- 

, lar cloud interfaces. In addition, S 106 IR drives an ionized wind 

■ with a velocity of -200 km s ' (ISimon & Fischeii Il982h that 
' is responsible for the hourglass shape and high velocity wings 

seen in optically thick molecular line emission close to the star. 
Here, we focus on the immediate surroundings of SI 06 IR, rel- 
evant for the context of the SOFIA observations . More details 
and references a re given in [^ odapp & SchneideH (l2008l) and by 
ISchneider et aP (120021 l2003l i2007i) . 

The two lobes of ionized gas seen in the optical and ra- 
dio continuum are separated by a dark lane that is very promi- 
nent i n near- to mid-IR imaging (ISmith et aI.L 1200 it lOasa et aU 
|2006|) . (S ub)Millimeter dust and optically thin molecula r line 
emission (iRicher et al.lll993HSchneider et aI.Ll200ll2003 l) show 
two emission peaks ~15" east and west of S106 IR in the dark 
lane, while the bulk of the molecular gas is located further to the 
east and west of the lobes. Two clusters of H2O masers were 



dete cted at the we s tern peak (S106 FI R, IStutzki et alL 119821: 
Fur uva etaU 1 19991) . [Richer et al.l (Il993h interpret this source as 
a Class young stellar object. 

While the dark lane w as initial l y inte r preted as a s mooth , 
large-scale di sk (Ballv & Scovillel Il982t iLittle etalL Il995h . 
iBarsonvetai] (IT989) and Richer et all (1 19931) concluded that 
the molecular line and dust continuum emission arises from a 
clumpy molecular cloud, possibly the dense remnants of a large 
disk or torus disrupted by the central star. A high degree of 
dumpiness was already inferred early on th rough observations 
of ammonia (I Stutzki & Winnewisseil Il985l) . Whether the dark 
lane is just a sha dow cast by a sma l l, edge-on ( accret ion) disk 
close to the star (' Ballv et aI.L Il983t iNoeletaLl l2005h . largely 
confining the ionizing radiation to the optical/radio lobes, or in 
addition due to dust extinction in high column density gas, is a 
matter of ongoing debate. A possible detection of the disk was 
reported by Hoare & Muxlow ( 1996) and Gibb & Hoare (2007'i). 

Observations in mid- to high-J CO lines, specifically tracing 
PDR gas close to SI 06 IR, indicate the presence of warm ( T>200 
K) and dense (n>10^ cm"^) gas ( Harris et al., 1987: Graf et al 
1993; Richer etal., 1993; Little et all 119951: ISchneider et al. 



i2003i) . iSchneider et al.i (i2003i) observed a larger area (few ar- 
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cmin) in [Cii] 158 yum using the Kuiper Airborne Observatory 
(KAO), and submm [Ci] and CO lines (KOSMA 3 m), reveal- 
ing the whole spatial extent of the PDR region. Higher angular 
resolution [O i] observations towards the center provide a more 
detailed view of the higher density PDR gas. The [C ii] and [O i] 
KAO data, however, were not velocity resolved. 



2. Observations 

The [Cii] atomic fine structure line at 1.900537 THz (157.737 
yum) and the CO 11^10 rotational line at 1.267 THz were ob- 
served with the heterodyne receiver GREAlQ on board SOFIA 
during two flights on April 13 and July 22, 2011, from 
Palmdale/California. Total power on-the-fly maps with a scan- 
ning speed of 8"/s, 1 s dump time, and a typical size of 2'x2' 
were performed. Positional offsets refer to the star S106 IR 
(R.A.,Dec.)(J2000)= (20''27"'26^74,37°22'47".9). The region 
centered on SI 06 IR has three coverages while the western 
extension has only one to two coverages, resulting in a non- 
uniform noise distribution. Blank sky subtraction was done to- 
wards a position off'set by (8',0'), well outs ide any emissio n 
seen in the KAO [C ii] map presented by Schnei der et al.l (l2003h . 
Instrument alignment and telescope efficiencies, antenna tem- 
perature and atmospheri c calibration, as well as the spectrome- 
ters us ed are described in lHevminck et al.l(l2012l) and lGuan et al.l 
(I2OI 2I). Here, we only show data from one of the Fast Fourier 
Transform Spectrometers (AFFTS), the other spectrometers giv- 
ing redundant information. All line intensities are reported as 
main beam temperatures scaled with main-beam efficiencies of 
0.51 and 0.54 for [C n] and CO, and a forward efliciency of 95%. 
The mean r.m.s. noise temperatures per 0.5 km s"' velocity bin 
for the central region (western extension) are 2.2 (3.7) K for [C 11] 
(16" beam) and 1.8 (2.4) K for CO (19.6" beam). The absolute 
calibration uncertainty is estimated to be ~10%. The IRAM 30 
m observations used for comparison w ith the new SOFIA data 
are described in lSchneider et al.l (12002). 

The 350 /xm continuum observations were obtained in 
March, 2004, with the SHARC-II bolometer camera at the 
CSCq Our scanning strategy was to modulate the telescope 
pointing with a non-connecting Lissajous pattern within the 
limits of reasonable telescope acceleration (typical periods of 
20 s). Calibration was frequently verified on planets and moons. 
Pointing was checked about every hour on evolved stars and 
blazars. Intensities are given in Jy/beam with a beam size of 9". 



3. Results and discussion 

3.1. Global morphology and the dark lane 

Figure [T| shows the integrated intensity distribution of [C 11] 
158 jum and CO 11^10 emission as observed with SOFIA. The 
emission morphology in both lines shows a pronounced peak 
centered on SI 06 IR and extended emission towards the south, 
north-west, and north-east, creating the heart-shaped emission 
already seen in the KAO data. The CO emission is more con- 
fined to the central region, offset from the emission peaks seen 



' German Receiver for Astronomy at Terahertz. GREAT is a devel- 
opment by the MPI fiir Radioastronomie and the KOSMA/Universitat 
zu Koln, in cooperation with the MPI fiir Sonnensystemforschung and 
the DLR Institut fiir Planetenforschung. 

^ The Caltech Submilhmeter Observatory is operated by the 
Cahfomia Institute of Technology under cooperative agreement with 
the National Science Foundation (AST-0838261). 



in [Cii], and better matches the [Oi] emission from the dense 
PDR gas JSchneider et all l2003l) . The SHARC-II 350 fim data 
show SI 06 FIR as a bright, unresolved source west of SI 06 IR, 
while east of the star, the emission is elongated to the north-east 
with some fainter extension towards the south-east resembling a 
dusty cone. Close to SI 06 IR, the integrated [Cii] intensity has 
two lobes perpendicular to the submm continuum emission. 

Near-IR and cm continuum show the northern and southern 
lobe divided by the dark cone in the east. Near SI 06 IR, the two 
[Cii] peaks closely follow the cm continuum, suggesting that 
part of the [C 11] emission is coming from the H 11 region, while 
the more extended [C 11] is unrelated to the cm continuum. The 
SUBARU image does not trace the full extent of the H 11 region 
just north- and south-east of S106 IR. Here, there is clearly ra- 
dio continuum emission from behind the dark cone in the fore- 
ground. The submm continuum emission perfectly matches the 
dark lane as well as the fainter near-IR dark regions just south of 
it. These observations reveal that the dark lane and cone are not 
just a shadow of the small disk around S 106 IR, but indeed high 
column density, warm gas at the edge of the molecular cloud 
(traced by the -2 km s ' channel of '^CO 2^1 in Fig. lA.lb . lo- 
cated in front of the H 11 region. The shadow of the small central 
disk protects the dark lane from ionizing UV radiation and likely 
is responsible for the fact that the lane survived for a significant 
time of the nebula's evolution after the formation of the central 
star. We will now use our new SOFIA and complementary IRAM 
30 m spectral line data to study the kinematics of the warm gas 
in the context of the nebula structure and geometry. For this pur- 
pose, we include various channel map overlays in the Appendix 
that reveal how complex the emission in the different velocity 
components and for the different tracers is. We, therefore, fo- 
cus this letter on the analysis and discussion of the morphology 
and kinematics of S106 and defer a more quantitative analysis, 
including PDR modelling, to a later paper. 

3.2. Velocity structure 

Channel maps of [Cii] and CO 11— >10 emission (top panel of 
Fig. lA. Il l reveal the different velocity components of warm and 
dense gas in S106. The emission can be divided into three ma- 
jor ranges for [Cii]: (1) localized blue emission close to SI 06 
IR (v=-6 to -3 km s"'), (2) extended emission with the charac- 
teristic heart-shape and an extension to the south at R.A. offset 
-70" associated with the velocity of the bulk of molecular gas 
(v=-3 to 1 km s ', see low-J CO spectra in Fig. |2]and chan- 
nel maps in Fig. lA.ll ). (3) and faint red emission spatially coin- 
ciding with the southern Hii region lobe (v=0 to 4 km s"' and 
higher). CO 1 1 —> 10 is offset from the [C 11] emission in all chan- 
nel maps, with [C 11] often wrapping around the CO. These differ- 
ences are due to the different excitation conditions. CO 11^10 
is prominent closer to the star as it traces the warmer, dense 
clumps mainly associated with the dark cone close to S106 IR 
(the J=l 1 level is 365 K above ground and the critical density is 
>10^ cm"-*), while [Cii] traces PDR or shocked material at the 
interfaces of the cloud and the H 11 region. 

The finding that the northern lobe is dominated by blue 
and the southern by red [Cii] emission is surprising as the 
southern lobe is tilted towards the observer while the north- 
ern lobe is pointing away (e.g.. ISolf & Carsentvl (119821) 1. For a 
proper explanation, we need to consider what is known about 
the geometry and dynamics of the nebula. The lobes of the 
hourglass nebula expa nd both radially and along the main axis 
dHodapp & Schneideit r2008) driven by S106 IR. The wind and 
the expanding lobes are expected to sweep up material where 



R. Simon et al.: SOFIA observations of S 106: Dynamics of the warm gas 




50 
20''27"'30 



-50 -100 

Offset (arcsec) 

20° 



-150 




50 
20''27"'30 



-50 -100 

Offset (arcsec) 

20° 



-50 

Offset (arcsec) 



-100 



-150 



-50 

Offset (arcsec) 



-100 



-150 




-150 



Fig. 1. Top: Integrated intensity (vlsr -20 to 20 km s"') of [C ii] (color) and CO 11— »10 (contours from 15 to 215 in steps of 15 K 
km s '). The position of S 106 IR is indicated by a star and the SOFIA beam size for [C ii] is drawn in the bottom left corner. On the 
right, contours of [C ii] integrated intensity are shown on the SHARC-II 350 jum image (intensities in Jy/beam). Bottom: Overlays 
of [C ii] on 1.47 GHz radio continuum (left, VLA archival data) and SHARC-II 350 /zm contours (right) on the SUBARU near-IR 
image. 



they hit molecular gas, except for the front wall of the south- 
ern lobe which has been partly eroded by the star (the back side 
of the southern lobe is visible in the optical and visual extinc- 
tion has been shown t o be significantly higher towards S106 IR 
and the northern lobe (iFelU et al.Lll984l) V We thus interpret the 
higher velocity blue and red [Cii] emission (v<-4 and v>l km 
s"') as arising from swept-up material at the front and back sides 
of the expanding, wind driven hourglass. In this picture, one ex- 
pects net blue shifted emission (with respect to the bulk molec- 
ular emission at — 1 km s'') from those parts of the hourglass 
close to S 106 IR that are curved towards the observer (front side) 
and red shifted emission from those curved away (back side). 
Indeed, blue shifted CO and [C ii] emission is particularly promi- 
nent in the northern and southern lobe close to SI 06 IR where 
the dark cone converges towards S106 IR in the foreground (-4 
to -2 km s"' velocity channels). 

As most of the front wall along the main axis of the southern 
lobe has been eroded, we do not expect extreme blue shifted 
emission there. The backside of the southern lobe should also 
show blue emission because of its tilt towards the observer. We 
suggest that this emission is shifted towards the red due to the 
radial expansion of the lobe. Red shifted emission is observed 
only towards the southern lobe very close to SI 06 IR and quite 
faint further along its axis, which may indicate that the back side 
of the northern lobe has also been eroded. 

For emission around the velocity of the bulk of the molec- 
ular gas, [Cii] traces the side walls of the lobes (i.e., where the 



radial component of the velocity is smallest). Comparing [C ii] 
and CO 1 1^10 emission to '^CO 2^1 (Fig. lAlT) allows to bet- 
ter characterize their relation to the bulk of the molecular cloud 
that is well traced in the low-J CO line. [C ii] (and to a certain de- 
gree CO 1 1— >10) in the -3 to +1 km s ' range nicely traces the 
surfaces of the molecular cloud. We therefore attribute the [C ii] 
emission in this velocity range as arising from PDRs in the cav- 
ity walls of the ionized lobes, directly correlated with the bulk 
emission of the molecular cloud. Further support for the above 
scenario comes from the overlays of [C ii] velocity channels on 
continuum data from the VLA (emission from the H ii region) 
near-IR from SUBARU (emission of hot dust), and SHARC-II 
(high columns of warm dust) in Fig. IA.2I The comparison with 
the 350 jum data is particularly interesting as it shows that [C ii] 
is confined to or even funneled into the lobes by the high col- 
umn density, warm dust seen in the submm continuum (the dark 
cone and S106 FIR). We even observe the dark lane sandwiched 
between [Cii] and '^CO 2— >1 emission, showing that there is a 
layering of warm ionized gas, warm dust, and colder molecular 
cloud material at the surface of the cloud. The strongest [C ii] 
emission originates from the edge of S106 FIR which is prob- 
ably a signature of ablation or evaporation of this dense clump 
due to radiation and/or wind from the star. 

Assuming optically thin, thermalized emission and beam fill- 
ing of unity, we estimate lower limits to the beam averaged [C ii] 
colu mn density for the em ission in the channel maps follow- 
ing ([Crawford et al.L Il985h in the high density, high temperature 
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Fig. 2. [Cii] and CO 11->10 spectra from SOFIA together with 
low- J CO lines from IRAM ( Schneider et al. , 2003 ) at a common 
angular resolution of 25". Offsets are given in arcsec from S106 
IR. From top to bottom: S106 IR, north-east lobe, and molecular 
cloud. 



limit. For the [C ii] wave length of 158 //m, an integrated inten- 
sity of 1 K km s"' translates into 7.041 x 10"^ erg s"' cm"^ sr"' 
and a column density of 4.6 x 10''' cm"^. The total [C ii] column 
density towards the emission peak is 2.0 x 10'** cm"^, while the 
values for the channel maps in 2 km s ' bins span a small range 
from 1.4 to 2.4 xlO'^ cm"^, the maximum being in the -A km 
s ' channel. The [C ii] column density in the red velocity bins is 
of the order of 7.3 x lO"" cm I 



3.3. Spectral signatures of high velocity gas 

The above findings are also reflected in the spectra displayed in 
Fig. 12] for three selected positions. The spectra have a common 
angular resolution of 25" and show line profiles towards SI 06 
IR, the north-eastern lobe (cavity wall), and the molecular cloud. 
[Cii] emission peaks at velocities of — 4 km s"', blue shifted 
from the peak velocity of the low-J CO lines at -1 km s ', and 
is generally very broad with emission covering the range from 
-14 to -1-12 km s"'. This implies that [Cii] traces very well the 
dynamics of warm gas in the PDR-layers, shocks, or the ionized 
phase, which most likely all have slightly different velocities and 
are subject to enhanced turbulence, quite diff'erent from the bulk 
of the molecular cloud. 

The CO 1 1 — > 10 line peaks near S 106 IR at velocities around 
-4 km s ' and shows broad lines due to the stellar wind hitting 
the molecular gas, thus tracing the warmest part of the molec- 



ular cloud edge in the foreground closest to the star where the 
wind speed is higher or even shocks can play a role. At this 
position, we sample emission from both the [C ii] and the CO 
11^10 peak. Towards the north-eastern lobe, CO 11— >10 be- 
comes much weaker and the peak velocity shifts towards that 
of the lower-J lines, likely reflecting that temperatures and den- 
sities are too low for sufficient excitation in the more quiescent 
gas outside the H ii region. Low-J CO lines with high to moderate 
optical depth also show blue wing emission, but not as extended 
to the far blue side as observed for [Cii] and CO 11— >10. This 
confirms that the latter lines trace gas in the foreground closer to 
the star and detached from the bulk of the colder molecular gas. 
High velocity red emission above -1-6 km s ' is only prominent 
in [C ii] and not seen in any of the other tracers. 

4. Summary 

Our new spectrally resolved observations of [C ii] and CO 
1 1— >10 reveal very complex morphology and kinematics of the 
warm gas in SI 06 with multiple velocity components at diff'er- 
ent locations relative to the exciting star, the dark lane, which we 
show to be composed of warm, high column density gas, and the 
bulk of the colder molecular gas. The [C ii] and high-J CO emit- 
ting material is highly dynamic with broad, non-Gaussian wings 
in particular in [C ii] showing no counterpart in any of the other 
observed tracers for the highest red and blue shifted velocities. 

Only the high spectral and angular resolution provided by 
GREAT on SOFIA make it possible to disentangle emission 
components arising from the PDR surfaces at the inner cavity 
walls of the two H ii region lobes, from swept-up gas due to the 
stellar wind of SI 06 IR and/or the expansion of the lobes, and 
from the H ii region itself. To quantify how much of the emission 
is contributed by shocks, H ii region, or PDR gas requires a more 
detailed analysis involving modelling, which will be addressed 
in a forthcoming paper. 
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Appendix A: Channel maps 




Fig.A.l. Channel maps of CO 11— >10 contours on [Cii] (top) 
and [Cii] and CO 11^10 contours on '^CO 2^1 (middle and 
bottom). 
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Fig. A.2. Top: channel maps of CII contours on a VLA 1.47 cm image. Middle: channel maps of CII contours on SUBARU near-IR. 
Bottom: channel maps of CII contours on SHARC-II submm continuum. 



